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A ciosed form solution is derived for the biphoton k-vector spectrum for an arbi- 
trary pump spatial mode. The resulting mode coefficients for the pump input that 
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resulting biphoton flux. 
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Oh- I. INTRODUCTION 



The effect of higher order or multiple pump spatial modes on the shape of the biphoton 
wavepacket resulting from spontaneous parameteric downconversion (SPDC) is of interest for 

generation of arbitrary single photon 



CN ! several reasons, including the determinsitic (heralded 
> 



> n n 

states [1[ 



entanglement of large numbers of photons [2|, super-resolving phase measurements 

and tl 

4- 



3|, and the entanglement of differing angular momentum states of the electromagnetic field 



Previous theoretical and experimental |7|, |8[ investigations into the spatial distri- 

bution of the quantum noise, in SPDC were rendered incomplete due to the inability to 



determine the biphoton spectra of higher order TEM modes. From a theoretical viewpoint, 
an experimentally testable form for the biphoton spectra of an arbitrary spatially shaped 
pump pulse has not been derived. Experimentally, the complete impact of differing experi- 
mental parameters on the corresponding parameters of the biphoton wavepacket is desired. 
We propose to solve these problems by choosing a different basis with a defined Fourier 
Transform and solving for a complete, closed form solution with Siegman's Elegant Gauss 
Hermite Modes, while only assuming that the width of the crystal when compared to the 
Guassian spot size is large such that the transverse Fourier Transform may be analyzed. 

In short, the "Elegant Gauss Hermite Modes," also known as the "Complex Gauss Her- 
mite Modes" (which we will use interchangeably throughout) have a defined, derivable 
Fourier Transform. This is due to the fact that they may be written as total derivatives 
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acting on the paraxial diffraction kernel. We will use this alternative mode expansion to 
convert the interaction Hamiltonian of a single nonlinear crystal from an integral over the 
crystal volume to an integral over the momenta of the signal and idler photons, thus de- 
riving the shape of the biphoton spectral density for an arbitrary superposition of Complex 
Gauss Hermite Modes. The form of the biphoton spectral density for arbitrary spatial pump 
distributions has not been derived to the author's knowledge. It is thus impossible to tailor 
a pump beam to optimize the biphoton distribution along anything but the optic axis of 
the pump laser without this solution. Although we will not do so here, in the far field, the 
biphoton spatial frequencies may be interpreted as the angle that the signal or idler photon 
makes with the pump beam optic axis; we may thus convert the integral over a correspond- 
ing detector area into an equivalent integral over the small region of spatial frequencies that 
are sampled by a detector with limited area. 

We will start by familiarizing the reader with the notation to follow. The notation is 
included in a helpful table. 

To simplify the notably variable laden integrals to follow, we will use £ = 1 + if-, 
following Boyd 9], and p 2 = x 2 + y 2 throughout. uo(t) is a normalization constant from 
Seigmann that will later represent the time variation of the amplitude (it will remain 
classical, as the pump field is treated classically; technically, it is "quantized" and we 
assume a coherent single mode, albeit Elegant Gauss-Hermite, and replace the raising and 
lowering operators with their eigenvalues) 10] . Sums are explicitly expressed, r_i_ is the 
component of r that is perpendicular to the optical axis (by convention, we define positive 
"z" to be the optical axis of the pump laser). Also, v + = v s + z/j, or the sum of the signal 
and idler fields that will be defined later. Any subscript "+" is to represent the sum of 
the signal and idler variables of that particular variables. We will also use e s to represent 
the two dimensional polarization vector of the electromagnetic field. We point out that 
choosing £ = 1 + ij- rather than £ = 1 — %— is equivalent to choosing the form k ■ r — ut 
for the phase of the electromagnetic field, rather than its negative. 
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An elegant Gauss-Hermite mode is denned as 



10|: 



u, 



U - 
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d n d n 



(1) 



10( to the form stated above is contained 



£ dxdy 
The derivation from Siegman's original form in 
in the Appendix. 

An elegant Gauss-Hermite mode of order n,m is biorthogonal with the function Vw : [ll3]: 

dxdy U nm 1p nm $nmfnm (^) (2) 

ip nm = H n [ —x ) H m [ —y ) (3) 



Xz r ^* 



Although the Hermite polynomials here are complex, and in fact the solution to a non- 
Hermitian differential equation, we may reasonably expect that a physicaly generated laser 
beam (such as might be expected for the pump beam) will show a distribution that asymp- 
totically converges to zero as the transverse area tends to infinity. This is equivalent to 
assuming that the total energy of the pump beam is finite. If this convergence is Gaussian, 
as is common, and we set the origin of the Complex Gauss Hermite Modes to the same origin 
as the pump beam, we can see that we may assume that the resulting polynomial expression 
in front of the Gaussian phase factor will be insufficient to disrupt the aforementioned con- 
vergence. In short, we can assume that, for a physical laser beam generated in a cavity (as 
is common for high power Q-switched lasers that are commonly used as the pump lasers in 
these experiments), the resulting spatial profile will be a superposition of "regular" Gauss 
Hermite modes, and the multiplication of two real Hermite polynomials is not sufficiently 
divergent to overcome the gaussian damping factor. This is possible because we choose the 
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origin of the Elegant Gauss Hermite Modes such that the overlap integral is performed at 
that origin, the only position where the modes must be real. The factor £ will be 1 at that 
origin. 

The fnm(z) is a constant when the origin of the elegant modes corresponds to the origin 



of the physical electric field 10]. 



f A-w ) n (-w ) m d n d m 
4 dxdy 

We may write the electric field in its entirety as, with 1 as the polarization index: 

E p {r, t) = e 2 ™^ Yl u nmi(t)cnmei + C.C. (5) 

This form of the electric field satisfies the paraxial wave equation. This can be seen by 
direct application of the paraxial wave equation. The paraxial wave equation is: 

^(r) = «(r)e ifc *Vlu(r) + 2 ^^p = (6) 

We now find the Fourier Transform of the elegant Gauss-Hermite modes; we list all 
relevant Fourier Transforms that are required to analyze the modes in this manner: 

oo 

3„ f(^\—2mu-T 



FT(f(r)) = F(u a , u v , v z ) = jjj d A x f{v)e~^ (7) 

r poo 

FT ± (f(r)) = F(u x , u y , z)= dxdy /(r) e - 2 — r (8) 

J J —CO 

FT± ^ ndmf lx^ {v,h) ) = ^(^n^ y rF(u x , u v ) (9) 

FTj_ ^ e -7r ( (:c/a)2+{s//b)2 ) \ = abe~^ {aUx)2+{bUv)2 ^ (10) 

(11) 

A brief note of importance here: first, ()9]) does not have the scaling factor that is normally 
present in the spatial frequencies; this is expected for two reasons: first, an analysis of the 
units associated with the Fourier Transform demonstrates that if the scaling factors were 
present, the result would be dimensionally incorrect. Second, the scaling factor cancels in 
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the following manner, which we will demonstrate in one dimension: 

x' = x/a (12) 
dx' = dx/a (13) 

dx' dx ^ ^ 

f d n \ C d n 

FT {dx^ fix/a) J = J dx dx^ f ^ a ^ X ^ 

(d n \ f d n 

^/(*/<>)J = a~" J (adx')— /(x')e i2 "C-»' (16) 

( d n \ 

FT l—f( x /a)\ = a{i2nv x ) n F(av x ) (18) 

This fact will allow us to completely derive the form of the interaction Hamiltonian in 
three dimensions, resulting in a derivable phasematching function (else the Fourier Transform 
in z would not have a closed form). 

Relating these to the Fourier Transform of the Complex Gauss Hermite modes, we see 

„,„, - I r((x/a) 2 + ( H /b) 2 ) 

that we have a function of the form e dx „ dyn , where a and b are, for a circular 



Complex Gauss Hermite mode, equal. They are both (for a circular case) equal to i/Az r £, 
and the multiplication of the two (ab) will give Az r £, which will, in combination with the 



, ^CLll^^l bLl^ Ull Lll^: U.^llVJlllllldUWl . llllO Will CL11WVV UO I 

longitudinal integral, which is merely a plane wave. Noting that Xz r is ttwq, or the transverse 



prefactor ( w °£ , cancel the £ on the denominator. This will allow us to evaluate a simple 



area of the beam, and collecting like terms, the Fourier transform is: 

FT(u nm ) = u (t)A ± (-2mw u x ) n (-2mw v y ) m e ~^(^) (19) 

Aj_ = ttWq is the transverse area of the beam (actually, this is the e~ 2 area of the beam, or 
approximately 95% of it). These results are expressed in a form sufficient for generalization 
to elliptical Elegant Gauss-Hermite modes (simply, the spot sizes in two dimensions need 
not be the same, and are merely scaling factors for the spatial frequencies in the appropriate 
dimension.). 

We will use ffTTl) in the derivation of the interaction Hamiltonian. 
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II. DERIVATION OF THE ELEGANT GAUSS HERMITE INTERACTION 

HAMILTONIAN 

The derivation of the interaction Hamiltonian will follow quite simply from the interaction 



Hamiltonian for three wave mixing [111 ]: 

i?i(r,t) = eo / <Pr X (ry>"''>E?(T,t)E' B >(T,t)Ef i (T,t) (20) 
Jv 

We will treat the time, and thus spectral, dependence of the field in a following section. 
For now, it is understood that we will only determine the spatial frequency dependence of 
the interaction. We will also treat the differing indices of refraction, inherently included 
in v = Vf s /n, where "fs" stands for "free space", for the calculation of the interaction 
Hamiltonian. 

For simplicity, we will treat the nonlinear crystal as cut for type-II (for example, o-eo) 
phase matching with a large, essentially infinite, surface area that is normal to the pump 
beam optical axis, but a finite length along the optic axis; the assumption of the infinite 
transverse area is our only assumption. Such an assumption is common in nonlinear optics 
9|. As such, we may assume that the raising and lowering operators of the signal and idler 
modes always commute. We may also drop the polarization vector, as the polarization of 
the signal and idler modes is now orthogonal, and set the overall phase of the pump mode 
to zero. 

We will also assume a coherent state for the pump field such that it may be treated 
classically. We also assume that the nonlinear coefficient does not vary over the crystal 
length for simplicity. 

We will write the fields as: 



(21) 



E s (r,0 = / ^V 7 /^ [ia s {v s )eJ 2 <^-^ + H.C] 

Ei(r,0 = / [ifli(^)e 8 e^- r -^ + H.C.] (22) 

(23) 

Where v is the spatial frequency, m is the index of polarization, and f is the frequency 
of the photon; as such, hf is the energy of the photon, where h is Planck's constant. The 
fields are quantized in free space, and the nonlinear crystal is treated as a perturbative 
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region of interaction. We note that we chose to operate in i/-space rather than k-space as 
normalization constants of the Fourier Transform are unity in z/-space, and this leads to a 
more elegant form of the interaction Hamiltonian. 

For clarity, we now drop the polarization index; we assume that the signal and idler fields 
are perpendicularly polarized (Type-II SPDC) and thus commute; it is possible to solve the 
problem with frames of vector spaces for identically polarized signal and idler fields with a 
frequency overlap, which would constitute an overcomplete basis. However, for the purposes 
of clarity, we demonstrate the (new) technique here without such complications. Explicitly, 



The quantum mechanical interaction Hami 



tonian of interest for SPDC is proportional 



to a\a\ and reduces to a Fourier Transform 11): 

n,m 

We now take the Fourier Transform over the transverse area to find: 
Hspdc = 2hxAu (t)^2c nm J dzSv^Vi < 



x (2iriw u +x ) n (27riw u +y ) m 



(25) 

Where we have made use of (TIT]) . The overall minus sign results from convention. 

Note here that we were careful to ensure that the complex phase of the pump mode, 
e +ik p z an( j g auss j an phase factor, was positive, in keeping with our convention of defining 
the phase of the lower operator of the field in the same way. This was a direct result of 
treating the state of the pump field semi-classically; the complete Hamiltonian would include 
a term proportional to a p a\a\, and the coherent state input would merely replace a p with 
the appropriately scaled amplitude of the pump. In keeping with our intuition, integration 
over all space would return a delta function where the sum of the signal and idler k-vectors 
along the optics axis must be identical to the initial pump mode optical axis k-vector. 

If we label the momentum mismatch (i — Az/^ ± — v +z ) = Au, we find evaluation over 
the longitudinal, and thus collinear with the optical axis, dimension, z, gives the expected 
phasematching function [9] $ = e l2lvAuAz sin 7 [^^ L - > , multiplied by the length of the interaction 
L, where Az is positive for a Gaussian mode with a focal mode that corresponds to the back 
of the crystal. We will call V the volume of the interaction, defined for the moment as AL, 



s 



where A is the transverse area of the beam, ttWq, and L is nominally the crystal length, 
although properly it is the length over which the interaction is coherent. 



Hspdc = -2¥%(t)Vc„ / d 3 u s d% { \ (26) 

III. TIME EVOLUTION OPERATOR AND THE BIPHOTON STATE 

In this section, we derive the time evolution operator to first order in the nonlinear 
coefficient %■ It is defined as 12]: 



U(t,t ) = ^(^j ±J*dt 1 £dt 2 ---£dt n T[H i (t 1 )H i (t 1 )---H i (t n )} (27) 
We insert the time dependence of the pump field into the amplitude tto(t); the time 



121 ] . With this substitution, we find 



dependence of the raising operators is well known 
HsPDcit) cx u (t)e i2 ^+ where /+ = f s (u s ) + /fa). 

We treat the nonlinear coefficient \ as small and expand the propagator to first order to 
find, noting the Fourier transform of Uo(t) is denoted as u(f): 



U(t,m) ~ 1 - ~ / dt'Hrtf) (28) 

'to 



C/(oo,-oo) Pa 2hxVu c nm / d 3 u s d 3 Uiu(f. 

n,m 



a\v s )a\vi)y/ ' f s (u s )fi(ui)(27riw u +x ) n 
x (27iiw v +y ) m $(Ais)e- nXzr ( u +* +u +y) 



(29) 

We may thus derive for the biphoton state vector, noting that the initial state is the 
vacuum state, and assuming the interaction is "turned on" for a time approaching infinity, 
the following state: 

|V>/(t = oo)) = U(oo, -oo) |0) (30) 

(2mw v +x ) n (2mw v +y ) m <S>(Av) 



\ipi(t = oo)) = 2hxVu c nm J d 3 Vgd 3 ViU(h 



Wi) \ u s) 



(31) 

With rare mathematical exceptions that might be separable, the multimode state is thus 
obviously entangled. We now proceed to derive a similar Hamiltonian, but include the sum 
over multiple polarizations and multiple crystals of the same length pumped by the same 
beam. 
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IV. GENERALIZED HAMILTONIAN FOR SPONTANEOUS 

DOWNCONVERSION 

If we wish to include polarization, we simply insert the polarization index into both \ 
and the polarization of the fields. We then find: 



Hspdc = -2/i ^2 c n m X° qr £po£* sq £ir Vu o J d 3 is s d 3 Vi 



n,m,o,q,r 



a \q{ v s)a\ r {vi) V fs(v s )fi(vi) (27riw u +x ) 

(32) 

Where e po , e* q , and e* r are the polarizations of the pump, signal, and idler respectively, 
summed over the indices o, q, and r. If we have multiple crystals, of the same length, 
being pumped by the same beam, we can write the spatial dependence of the interaction 
Hamiltonian as the "regular" interaction Hamiltonian as convolution with a sum of delta 
functions. We point out that this results merely in the sum of the individual interaction 
Hamiltonians, and we are considering a simplified case with some interesting results. 



V. MAXIMIZATION OF THE BIPHOTON SPECTRUM FOR A DETERMINED 

SOLID ANGLE 

In this section, we maximize the biphoton spectrum for a desired set of signal and idler k- 
vectors, and find the corresponding optimal expansion in the Elegant Gauss- Hermite Modes. 
We will do this by invoking a variational principle and determining the structure of the 
pump spatial mode that will produce the largest probability of measuring a biphoton for an 
arbitrary set of signal and idler k- vectors. As a result, we will also show that the only way to 
produce signal and idler k-vectors along the optic axis is to use the zeroth order Transverse 
Electromagnetic mode (we note that this is identical in the case of the CGH and regular 
GH modes). 

We start with our interaction Hamiltonian, assuming an initial vacuum mode for the 
signal and idler photon numbers: 

aL(vs)a\ r (iSi)^ f s (v s )fi(vi)(2niw is +x ) n 



Hspdc = -2h ^ c nm x° qr e po e* sq e* r V J d z v s d?Vi 



n,m,o,q,r 



x (2mw u +y ) m ^(Au)e- 7rXzr (' / ^ +u +y) 

(33) 
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The integration will be performed over the detector volume, which we will assume is in 
the far field. We thus wish to maximize the probability of the biphoton arriving at the 
detector. We see that, assuming the interaction Hamiltonian is small and that the resulting 
measurement will be a two-photodiode correlation in the far field, which we may express as 
a finite integral over the k-vectors that the detectors are sensitive to, and approximate the 
resulting integral as slowly varying in the far field, and proportional to the functional form 
of the interaction Hamiltonian multiplied by the accompanying sensitive k-space volume of 
the detector (which will be labelled A(V S ) 3 and A(Vj) 3 ), the resulting term, state vector of 
interest will be proportional to: 



m=H SP Dc\0) (34) 
Hspdc = -2h c n m X oqr ^poe* sq e* r Vu J Sv s d 



n,m,o,q,r 




a lq ( V s ) 4r ( V i ) \J f 8 ( V s ) f { ( V { ) ( 2 TXIWq V +x ) 

x (2mw u +y ) m ^(Au)e- 7rXzr (' / ^ +u +y) 

(35) 



Hspdc ~2h ^ c nmX oqr e po e* sq 6* r Vu A(v s ) 3 A\ 



n,m,o,q,r 




( v s )a\Xvi)\J f s {v s )fi{vi){2niwov +x ) n 

(36) 

Where the v s and terms are now a single value. 

Any measurement that is to maximize the biphoton probability must maximize the cor- 
responding bilinear term in Hspdc- If we dub our measurement operator "M", we find 
that: 

M oc J d 3 V(ip | a\a\a s ai \ ip) (37) 
M oc Yl < m c ij(-^ok +x ) m (-tw k +y ) n (iw k +x )\iw k +y y)) (38) 

nm ij 

We thus wish to maximize the last term. As we are working in the paraxial approximation, 
we may assume that Wok +x <<< 1, and a similar assumption may be made in the y direction. 
The maximization problem we must solve is also subject to one constraint, namely that the 
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sum of all the pump spatial mode coefficients must be 1. Explicitly: 



I Coo 1 2 + y ] \c n m\ 2 — 1 

CooCqo = 1 - Yl 



n—l.m—1 



(39) 

(40) 

(41) 

(42) 

(43) 
(44) 

We now minimize M, treating the mode coeffi cents and their complex conjugates as inde- 
pendent variables; this is equivalent to treating the real and imaginary parts as independent, 
as we may rewrite the mode coefficients as linear superpositions of the real and imaginary 
parts. Explicityly, we institute: = 0. 

We now proceed to find the critical points of M , noting that the sums over n,m, i, and j 
run from 1 to the maximum: 



c oo — 


1 Snm \ C nm\ 




coo 




Br* 


1 — Snm \ C nm\ 


c oo 


dc 00 


r 2 
c oo 


c 00 


Br* 


c* 




dc nm 


coo 





M 
dM 



E 



mm] 



dc 



c oo + c nm(- iw ok+ x ) m (iw k +y ) n )) (c 00 + Cijiiwok^yiiwok+yY)) (45) 

(46) 



on 



mni 3 C °° 



We may set the phase of c 00 to zero without loss of generality. As a consequence, we 
see that the term c* m ((— iwok +x ) m (iw k +y ) n )) must be real; this implies that we must set 
the phase of the modes to be (n + m)(— vr/2); this phase makes all terms in the interaction 
Hamiltonian real, and we may solve a simpler minimization problem; instituting the equaltiy 



\n+m ~ 



we see that our problem now reduces to maximizing a real function 



squared, which is equivalent to maximizing the function: 

J^c nm ((w k +x ) m ){w kl y )) 



M 



(47) 
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We now merely have to solve for the case of optimizing the term in the brackets to find: 

(48) 
(49) 
(50) 

(51) 

(52) 
(53) 

4,(l + ^X 2m F 2n ) = 1 (54) 



w k +x 


= X 


w k +y 


= Y 


Cnn 
^uu 


= sartl — 7 c 2 

" / j nm 

nm 


w 


= ^~c nm X n Y m ^ 


dM 
dcij 


mn 

= X l Y j - Q./Coo 


c i.i 





y /—/nm 



With this solution, we now see that the optimal coefficients between the modes is: 

= (-iw k +x y(-iwok +y y 
° 13 ~ V 1 + Enm(^k +x ) 2m (w k +y )^ 

We have thus maximized the biphoton spectra for a given set of k- vectors for signal and 
idler. We would also point out that the result will hold irregardless of the limits of the 
integral in k-vector space, as the minimiziation was performed over the mode coefficients. 



VI. CONCLUSION 



We have, for the first time (to the author's knowledge), determined the effect of multiple 
pump spatial modes on the resulting spectra of biphotons that result from the second order 
nonlinear spontanous interaction ("downconversion"). We achieved this result by expanding 
the incoming pump field (which we treated as classical, or equivalently assumed was in a 
coherent state) in a basis of Siegman's Complex Gauss Hermite modes, which, as we have 
shown, have a defined Fourier Transform. This results in, besides the normal gaussian phase 
factor that accompanies the Fresnel approximation in optics, a polynomial expansion in 



13 



the sum of the signal and idler frequencies. The form of the resulting biphoton spectra for 
type-II downconversion (e-oe or o-oe) is shown in [3TJ 
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